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ABSTRACT 


The  effect  of  signal  and  background  shot  noise,  as  well  as  device  noise,  on  the  performance 
of  a  direct-detection  spatial  tracking  system  is  investigated  for  arbitrary'  detector  arrays 
assuming  linear  loop  operation.  The  performance  of  quadrant  detectors  as  a  function  of 
background  radiation  and  detector  radius  is  then  analyzed  and  compared  with  some 
performance  bounds.  APD-  and  PIN-based  tracking  systems  also  are  compared.  The  effects 
of  non-focal-plane  processing,  focus  error,  and  pupil  walk  on  tracking  performance  are 
investigated.  Experimental  results  of  a  2.4-kHz  two-axis  tracking  loop  operating  at  low 
signal  power  using  an  APD  quadrant  detector  are  presented. 
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DESIGN  CONSIDERATIONS  AND  EXPERIMENTAL  RESULTS 
FOR  DIRECT-DETECTION  SPATIAL  TRACKING  SYSTEMS 


L  INTRODUCTION 

The  spatial  tracking  system  for  optical  radar  and  communication  is  of  critical  importance  because 
optical  antennas  have  such  narrow  beamwidths  (~|irad).  Even  a  ver>'  good  optical  platform  may  have 
angular  disturbances  on  the  order  of  a  few  milliradians,  and  optical  systems  typically  require  angular 
stabilization  to  a  fraction  of  a  beamwidth.  Thus,  a  spatial  tracking  system  with  high  precision,  large 
bandwidth,  and  large  dynamic  range  is  necessary.  A  number  of  papers  have  been  written  on  the  design 
of  such  tracking  systems.'"^  Here,  we  discuss  some  additional  design  considerations  and  present  some 
experimental  results. 

A  simplified  block  diagram  of  a  typical  spatial  tracking  system  is  shown  in  Figure  1 .  The  tracking 
system  control  signal  is  generated  from  the  tracking  detector  array,  which  is  typically  located  in  a  focal 
plane.  If  the  angle  of  incidence  of  the  received  optical  field  is  accurately  tracked,  the  focal  spot  will 
remain  centered.  If  the  spot  is  not  centered,  error  signals  are  generated  and  fed  back  to  a  fast  steering 
mirror  (FSM)  to  correct  for  the  errors.  Typically,  to  achieve  high  bandwidths  the  FSM  has  a  small 
angular  range.  Therefore,  to  keep  the  FSM  within  its  operation  range,  the  low-frequency  large-amplitude 
disturbances  are  nulled  with  a  large  coarse  pointing  mirror  (CPM)  or  a  gimballed  telescope.  This  nested 
servo  design  can  provide  a  dynamic  range  of  more  than  10^:1,  low-frequency  disturbance  rejection  in 
excess  of  140  dB,  and  -3-dB  closed-loop  bandwidths  in  excess  of  2.4  kHz. 


Figure  1.  Tracking  system  structure. 
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In  this  report,  we  focus  on  the  high-bandwidlh  tracking  loop.  We  develop  bounds  and  estimates  of 
noise  performance;  compare  the  performance  of  two  commonly  used  detectors,  avalanche  pholodeleciors 
(APDs)  and  PIN  detectors;  discuss  the  effects  of  wavefront  error  and  non-focal-plane  processing;  and,  in 
the  section  on  experimental  results,  discuss  a  iracking-error-normalizaiion  and  loop-compensation  tech¬ 
nique  and  present  quantitative  data  on  loop  rejection  and  noise  performance. 

In  optical  radar  and  communication  systems  the  boresighl,  acquisition,  and  pointing  subsystems  are 
closely  integrated  with  the  spatial  tracking  subsystem.  However,  these  subsystems  will  not  be  discussed. 
It  is  assumed  that  acquisition  has  been  accomplished  and  that  the  incident  optical  field  has  a  flat  phase 
front. 

2.  NOISE  PERFORMANCE 

2.1.  General  Noise  Performance  Bounds 

As  long  as  the  tracking  loop  operates  in  its  linear  region,  it  can  be  represented  by  the  standard 
model  of  Figure  2.  represents  the  angle  discriminator  gain.  In  this  simplified  model,  the  total 
tracking  error  is  given  by 

(3^  =  a^+NEA^  (1) 


where  is  the  uncompensated  tracking  error,  and  NEA  (noise-equivalent  angle)  is  noise-induced  track¬ 
ing  error.  Rewriting  this  expression  in  terms  of  the  power  spectrum  of  the  angular  disturbance  S^(0.  the 
closed-loop  transfer  function  of  the  tracking  loop  H(0,  and  the  noise-equivalent  spectral  density  (NESD), 
yields  the  well-known  expression 

oo  oo 

(4  =  j  I  1  -  H(f)  P  SgCf)  df  +  j  I  H(f)  1^  df  (2) 


ADDITIVE 

WHITE 

NOISE 


Figure  2.  Simplified  one-dimensional  linear  loop  model. 
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The  NESD  is  equal  to  the  (single-sided)  spectral  height  of  the  additive  white  noise  divided  by  the  square 
of  the  angle  discriminator  gain.  The  additive  noise  has  been  assumed  to  be  spectrally  white.  Most  of  the 
noise  sources  such  as  signal,  background,  and  dark  current  shot  noise  as  w^ell  as  thermal  noise  are  well 
modeled  as  white  noise. ^  There  are  circuit  and  detector  noise  sources  that  give  rise  to  non-w'hitc  noise, 
but,  as  we  will  see  in  Section  4,  over  a  small  bandwidth  the  approximation  of  a  white-noise  source  yields 
good  experimental  results.  Rewriting  Equation  (2)  in  terms  of  the  loop  rejection  R(0  and  noise  equiva¬ 
lent  bandwidth  (NEB) 

oo 

=  j  I  R(f)  r  SgCf)  df  +  NESD  NEB  (3) 


Therefore,  NEA  for  each  angular  axis  is  equal  to  the  square  root  of  the  product  of  the  NESD  and  the 
NEB.  The  noise  bandwidth  of  the  tracking  loop  as  well  as  the  uncompensated  tracking  error  are  gov¬ 
erned  by  the  dynamics  of  the  loop  and  angular  disturbance.  In  this  section,  we  develop  bounds  on  noise- 
induced  tracking  error.  To  make  these  bounds  independent  of  loop  dynamics  and  valid  for  arbitrary' 
angular  disturbances,  we  focus  on  minimizing  the  NESD  for  an  arbitrary'  number  of  detectors  placed 
symmetrically  about  the  x  and  y  axes.  In  Section  4,  we  will  discuss  tracking  loop  rejection  and  loop 
dynamics. 

To  minimize  the  overall  tracking  error,  it  is  necessary  to  minimize  the  NESD.  The  loop  dynamics 
can  be  optimized  separately  once  the  angular  disturbance  is  modeled.  It  can  be  shown  that,  for  a  linear 
tracking  loop,  minimizing  the  NESD  is  necessary  to  achieve  a  minimum  mean-square  error.  Further¬ 
more,  as  the  number  of  photons  per  coherence  lime  of  the  angular  process  becomes  large,  the  Poisson 
statistics  governing  the  observation  equation  approach  Gaussian.  Minimizing  the  NESD  is  then  neces¬ 
sary  to  achieve  many  other  performance  bounds  such  as  the  Cramer-Rao  bound.  * '  It  should  be  noted  that 
for  certain  loop  designs  the  uncompensated  tracking  error  may  dominate  the  NEA.  Under  such  circum¬ 
stances,  optimizing  the  NESD  will  not  dramatically  reduce  the  overall  tracking  error  but  will  add  addi¬ 
tional  margin  to  the  tracking  system  signal  power  link  budget. 

Assuming  that  the  loop  is  operating  in  its  linear  region,  the  currents  out  of  each  detector  element 
can  be  linearized  about  the  azimuth  and  elevation  tracking  errors.  Expressions  for  the  noise  out  of  each 
detector  based  on  signal  shot  noise,  background  shot  noise,  dark  current,  front-end  amplifier  noise,  and 
excess  noise  factor  (in  the  case  of  detectors  with  gain)  can  be  derived.*^  Optimally  combining  all  the 
detector  outputs  yields  the  minimum  NESD,  which  is  given  by** 


NESD=  [ 
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=  4  J*  Ijj(x,y)  dxdy 
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p 

m 

= 

normalized  intensity  modulation  power 

P 

s 

= 

optical  signal  power 

■n 

detector  quantum  efficiency 

V 

= 

optical  frequency 

h 

= 

Planck’s  constant 

I,(x,y) 

= 

normalized  signal  intensity  distribution 

Vx,y) 

= 

background  intensity  distribution 

x  and  y 

represent  spatial  coordinates  at  the  detector 

Q 

integration  area  across  a  detector  element 

F(G) 

= 

excess  noise  factor 

^db 

= 

gain  dependent  dark  current 

^ds 

gain  independent  dark  current 

e 

= 

electron  charge 

N 

oc 

= 

total  single-sided  front-end  current  noise  ! 

G 

detector  gain. 

The  infinite  summation  is  over  one  quadrant,  and  the  subscripts  “ij”  identify  parameters  associated  with 
each  detector  element. 

In  many  tracking  systems,  the  use  of  intensity  modulation  is  warranted  because  it  eliminates  sensi¬ 
tivity  to  1/f  and  dc  disturbances  from  the  detectors,  amplifiers,  and  background  sources.  Since  the  mod¬ 
ulation  depth  cannot  be  >1,  must  be  <1.  Assuming  coherent  demodulation  and  a  modulation  depth 
of  m,  P  =  m^  for  square-wave  modulation,  =  (4/k)"  m^/2  for  square-wave  modulation  using  only  the 
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fundamental  harmonic  in  demodulation,  and  =  m^/2  for  sinusoidal  intensity  modulation.  For  ideal  C\V 
operation,  =  1  • 

In  addition  to  increased  modulation  power,  square-wave  modulation  in  some  cases  offers  an  addi¬ 
tional  performance  improvement  in  the  presence  of  background  radiation.  Such  cases  arise  when  the 
optical  frequency  is  modulated  by  the  same  parameter  that  is  varied  to  modulate  the  intensity  (i.e., 
semiconductor  lasers).  With  square-wave  modulation,  an  interference  filter  need  only  be  w  ide  enough  to 
encompass  the  optical  frequency  uncertainty  when  the  source  is  in  its  high-power  state  (~A);  whereas, 
with  sinusoidal  modulation,  the  interference-filter  bandwidth  must  be  wider  than  the  frequency  uncer¬ 
tainty  plus  the  modulation  bandwidth  (~10  A).  A  narrower  interference  filter  leads  to  reduced  back¬ 
ground  noise. 

Equation  (4)  can  be  optimized  for  a  wide  variety  of  parameters.  To  compare  our  results  w  ith  those 
of  others,  we  analyze  Equation  (4)  for  three  different  cases. 


2,2.  Signal  or  Background  Shot-Noise-Limited  Operation 


It  is  assumed  that  the  detector  dark  currents  and  front-end  amplifier  noise  are  equal  to  zero  (G  =  1, 
F  =  1 ).  Under  these  assumptions,  there  is  no  loss  in  assuming  that  each  detector  element  is  infinitesimal 
since  they  can  always  be  grouped  to  form  a  larger  detector  with  no  additional  noise.  In  the  limit  as  each 
Q..  approaches  zero.  Equation  (4)  becomes 


NESD  =  [  2 


hv 
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.  ^Is(x,y)  2 

^  dx  ^ 


0  0  Ij.(x,y)  + 


dx  dy  ] 


(8) 


Rewriting  Equation  (8)  in  terms  of  the  amplitude  of  the  signal  field  and  assuming  that  there  is  no  back¬ 
ground  radiation  yields  a  lower  bound  for  the  NESD,  which  is  given  by 


NESD  >  [  8  P. 


oooo 


(9) 


This  expression  is  identical  to  the  corresponding  bound  for  heterodyne  detection  except  for  a  factor  of 
2P^  (see  Reference  5).  That  is,  for  ideal  CW  operation,  direct  detection  is  superior  by  3  dB.  (This 
analysis  is  based  on  the  first  and  second  moments  of  the  observation.  As  pointed  out  in  Section  2.1,  the 
direct-detection  and  heterodyne-detection  have  similar  statistics  only  when  the  number  of  detected 
photons  per  loop  time-constant  becomes  large.) 

The  bounds  for  Airy-disk  and  Gaussian  focal-plane  intensity  distributions  are  given  by 

NESD>[P^^^^r'  Airydisk  (10) 

nv  2 
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NESD>[  P„^-L]''  Gaussian  (11) 

hv  ^ 

where  X  is  the  optical  wavelength,  d  is  the  aperture  diameter,  and  a  is  the  radius  at  which  the  intensity  is 
reduced  by  1/e  point.  (Note  that  a  is  related  to  the  radius  at  which  the  amplitude  is  down  by  1/e  in 
the  pupil-plane,  by  a  =  (1/V271)  (X/cd^). 

2.3.  Signal  or  Detector  Shot-Noise-Limited  Operation  (Gaussian  Field) 

It  is  assumed  that  there  is  a  Gaussian  intensity  distribution,  no  background  radiation,  unity  gain  de¬ 
tectors  (G  =  1,  F  =  1 ),  and  square  detector  elements  with  equal  area  and  distance  a  on  each  side.  The 
objective  is  to  find  the  NESD  as  a  function  of  a.  The  result  is  given  in  Figure  3,  where  X^  = 

X^  =  [(lye)  +  N^/2e“].  It  is  clear  that  there  is  no  merit  in  using  detector  arrays  larger  than  a  quadrant 
detector  except  when  the  device  noise  is  near  zero.  This  reflects  the  fact  that  spreading  the  signal  field 
over  many  detectors  to  increase  resolution  does  not  offset  the  increase  in  noise  from  the  added  detectors. 
Even  when  there  is  no  detector  noise,  the  improvement  of  an  infinite  detector  array  over  a  quadrant  de¬ 
tector  is  <2  dB.  As  we  will  see  in  Section  2.4,  for  an  Airy-disk  intensity  distribution  an  infinite  detector 
array  can  improve  performance  over  a  quadrant  detector  array  by  -4.8  dB.  Still,  in  many  applications  the 
increased  performance  afforded  by  higher  order  detector  arrays  is  not  worth  the  added  complexity. 

As  seen  in  Figure  3,  once  the  quantity  a/a  is  greater  than  ~2,  the  tracking  f>erformance  remains  rela¬ 
tively  constant  and  close  to  optimal.  It  is  of  interest  to  compare  the  tracking  results  with  those  from  esti¬ 
mating  the  centroid  of  the  intensity  profile.*^  Under  similar  assumptions,  the  optimal  a/a  ratio  for 
centroiding  was  found  to  be  between  2.1  and  3.5.  The  upper  limit  results  from  the  spot  size  becoming 
very  small;  it  can  be  lost  within  one  detector  element  unless  it  is  near  the  edge.  This  is  not  a  problem 


Figure  3.  NESD  as  a  function  of  ratio  of  detector  to  spot  size. 


6 


116455-3 


116455-4 


in  the  tracking  scenario,  since  the  spot  is  being  stabilized  at  the  intersection  of  four  detectors.  In  situa¬ 
tions  where  both  acquisition  and  tracking  are  to  be  performed  with  the  same  detector  array,  if  the  quan¬ 
tity  a/a  remains  within  the  intersection  of  the  optimal  tracking  range  and  the  optimal  acquisition  range, 
which  is  a/a  between  2.1  to  3.5,  near-optimal  tracking  and  acquisition  performance  can  be  achieved 
simultaneously. 


2.4.  Signal  or  Background  Shot-Noise-Limited  Operation  with  a  Quadrant  Detector 


Section  2.3  showed  that,  for  a  Gaussian  intensity  distribution,  a  quadrant  detector  was  essentially 
optimum.  We  will  now  find  the  optimum  radius  for  a  circular  quadrant  detector.  For  a  quadrant  detector, 
P^can  be  expressed  as  shown  in  Equation  (12),  where  r  is  the  radius  of  the  detector  and  \/d  represents 
the  diffraction-limited  beamwidth.  We  assume  that  the  background  radiation  is  from  an  extended  source 
with  a  spectral  radiance  of  N(\)  and  that  the  optical  bandwidth  is  W: 


Pk=4Ph  [ 

b  bo 


i‘ 


(12) 


Pbo  =  (frNa)?i^  w  (13) 

Assuming  an  Airy-disk  intensity  distribution  and  the  detector  dark  current  and  front-end  amplifier 
noise  equal  to  zero  (F  =  1 ),  the  NESD  as  a  function  of  detector  radius  can  be  found  (Figure  4).  For  the 
case  of  no  background  radiation,  the  optimum  detector  radius  is  about  1.22  X/d  (the  first  null  of  the  Air>' 
disk)  and  the  i:>erformance  is  4.8  dB  below  the  optimal  performance  given  by  Equation  (10).  As  the 
background  radiation  increases,  the  optimum  radius  slowly  decreases.  In  practice,  larger  FOV  detectors 
with  their  associated  larger  pull-in  range  are  often  used  to  minimize  the  requirements  of  the  spatial 
acquisition  system. 


Figure  4.  NESD  vs  quadrant  detector  radius. 
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2,5.  Comparison  of  APD-  and  PIN-Based  Tracking  Systems 


We  now  compare  the  performance  of  APD-  and  PIN-based  tracking  systems.  For  the  APD-based 
tracking  system,  the  gain  can  be  optimized  to  minimize  the  NESD.  For  some  APDs,  the  dependence  of 
excess-noise  factor  on  gain  is  given  by^'^^ 

F  =  (2-l)(  1  -k)  +  kG  (14) 

u 


where  k  is  the  effective  ionization  rate  ratio.  Assuming  small  k  and  large  G,  the  optimum  APD  gain  is 
given  by 


Go  =  [ 


I  N 

8(-^  +  il2£ 

"  2e2 


,1/3 


k(s2^. 

hv 


TIP.  I.. 

_ ^  4  ) 

hv  e  ^ 


(15) 


Figure  5  contains  the  comparison  of  each  tracking  system  to  an  “ideal”  detector  which  is  assumed 
to  have  the  unity  quantum  efficiency,  no  dark-current  noise,  and  no  front-end  circuit  noise.  The  parame¬ 
ters  assumed  in  the  comparison  are  given  in  Table  I.  Three  different  background  sources  are  investi¬ 
gated:  no  background,  strong  earth  background,  and  solar  background.  An  Airy-disk  focal-plane  inten¬ 
sity  distribution  and  a  quadrant  detector  with  a  radius  of  2.5  beamwidths  are  assumed.  The  APD  and 
front-end  noises  are  values  that  have  been  measured  in  our  laboratory.'^  The  comparison  is  made  in 
terms  of  the  ratio  of  the  power  required  of  each  tracking  system  to  the  “ideal”  detector  tracking  system  in 


SIGNAL  POWER  (photons/s) 


Figure  5.  APD  and  PIN  vs  ideal  detector  performance. 
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TABLE  1 

APD  and  PIN  Performance  Parameters 

X 

0.8  pm 

Detector  Radius 

2.5  X/d 

S,,  =  0.92 

S„  =  1.64/(X/d) 

N 

oc 

2.0  X  10'^^  A^/Hz 

Pb 

1 .83  pW  [N(X)  =  0.024  W/(m^A  sr)]  earth’'^'^^* 

140  nW  [N(X)=  1800  W/(m^A  sr)]  sun^'’’’^* 

w 

25  A 

APD  Parameters 

PIN  Parameters 

II 

O 

bo 

II 

O 

bo 

'db  =  0-01  pA 

o. 

CT 

II 

o 

I .  =1.0  nA 
as 

1  =  1.0  nA 

ds 

k  =  0.01 

k  =  0 
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*  Background  assumes  3-dB  loss  due  to  a  polarizer,  and  2-dB  loss  due  to  optical  attenuation. 

order  to  achieve  the  same  NESD.  As  expected,  for  background-limited  operation  (such  as  the  solar 
background)  or  for  signal  shot-noise-limited  operation  (>10'‘  photons/s),  the  APD  and  PIN  systems 
perform  identically.  However,  for  low  backgrounds  the  APD-based  system  can  offer  over  20  dB  of 
improved  performance.  The  plateaus  at  low  signal  powers  are  the  regions  in  which  the  PIN-  and 
APD-based  systems  are  background  shot-noise-limited.  Note  that  the  APD  performs  within  3  to  5  dB  of 
the  ideal  detector  unless  it  becomes  gain-dependent  dark-current-limited;  3  dB  of  this  loss  can  be  attrib¬ 
uted  directly  to  the  APD  excess-noise  factor  [Equation  (14)].  The  fact  that  the  additional  loss  is  small  is 
a  result  of  the  relatively  low  gain-dependent  dark  current  and  the  low  ionization  rate  ratio. 

3.  NON-FOCAL-PLANE  PROCESSING 

We  now  investigate  the  effect  of  non-focal-plane  processing,  focus  error,  and  pupil  walk  on  track¬ 
ing  performance,  assuming  that  the  incident  optical  field  is  Gaussian  and  that  an  infinite  area  quadrant 
detector  is  used.  Using  the  Fresnel  diffraction  equation  to  find  the  field  distribution  at  an  arbitrary  dis¬ 
tance  from  a  focusing  lens,  the  sum-channel  and  angle-discriminator  profiles  can  be  found  and  are  given 
by 

I  =  Ps  (’6) 
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Ax 


where 

(f )(  I  -7^) 

x  =  [i  +  [  ^  ^  n'”  *'*' 


-L  =  l  X 

Lo  f’Ro  (19) 

Z  =  sum-channel  profile, 

=  x-axis  difference-channel  profile, 
erf  (.)  =  error  function, 

0^  =  x-axis  angle  of  incidence  of  the  optical  field  at  the  focusing  lens, 

L  =  distance  from  the  focusing  lens  to  the  detector  array, 
f  =  focal  length, 

=  x-axis  lateral  offset  of  the  optical  beam  at  the  focusing  lens, 

=  radius  of  curvature  of  the  field  incident  onto  the  focusing  lens. 

We  now  investigate  the  sensitivity  of  the  NESD  to  the  placement  of  the  tracking  detector  array.  By 
assuming  small  tracking  errors  so  that  the  error  function  can  be  linearized,  minimizing  the  NESD  is 
equivalent  to  maximizing  the  slope  of  the  angle  discriminator  or  equivalently  minimizing  X.  Optimum 
performance  is  obtained  when  the  detector  array  is  placed  at  L^.  This  distance  is  equal  to  the  focal 
length  only  when  the  incident  optical  field  is  perfectly  collimated  If  L  =  +  €,  then 

The  first  term  on  the  right-hand  side  can  be  interpreted  as  the  diffraction-limited  spot  size,  the  second 
term  as  the  f-number  of  the  focusing  lens,  and  the  third  term  as  a  loss  factor.  For  example,  if  only  a 
1.5-dB  loss  in  tracking  performance  (X  =  is  acceptable  and  an  f#  =  1  beam  is  used,  the  tracking- 
detector  array  must  be  placed  to  within  one  diffraction-limited  spot  size  of  the  true  location  of  the  focal 
plane.  This  restriction  can  place  stringent  thermal  and  mechanical  tolerances  on  the  optical  system 
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design.  For  large  f-numbers,  the  linearization  used  to  obtain  Equation  (20)  ma\  be  invalid  and  Equa¬ 
tions  (18)  and  (19)  must  be  used  directly. 

We  now  investigate  the  sensitivity  of  the  NESD  to  focus  error,  assuming  the  tracking  detector  is 
placed  in  the  optimal  location  for  an  unaberrated  system  (L  =  f).  If  focus  errors  are  introduced  from 
somewhere  within  the  optical  system,  then  the  tracking  performance  will  be  degraded  as  a  result  of 
thermal  or  mechanical  distortion  or  a  design  residual.  If  N  is  the  number  of  waves  of  defocus,  then 


N  = 


1 

cr 


1  \n 

(x^-1 

4k 


(21) 


Again,  if  the  loss  in  tracking  performance  is  to  be  <1.5  dB,  the  focus  error  must  be  kept  to  <X/13. 

We  now  investigate  the  effect  of  the  lateral  offset  from  the  optic  axis.  As  seen  from  Equation  (17), 
if  the  detector  array  is  not  placed  at  the  focal  distance  the  offset  appears  as  a  bias  term  in  the  error  signal. 
If  the  bias  becomes  significant  compared  with  X/(2Ka)  X,  then  the  angle  discriminator  will  become 
nonsymmetrical  and  its  effective  slope  will  be  reduced.  The  dynamic  range  over  which  the  angle  di.s- 
crimination  remains  linear  can  be  chosen  by  adjusting  the  location  of  the  detector  L.  For  example,  if  a 
1-dB  compression  in  the  discriminator  linearity  is  acceptable  at  a  maximum  angle  then  the  detector 
array  should  be  placed  at 


L={ 


“  CJ  ^^0.875 
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However,  this  increased  dynamic  range  is  not  only  at  the  expense  of  noise  performance.  It  can  lead  to  a 
limited  angular  rejection  capability.  If  the  pupils  within  the  optical  system  are  not  tightly  controlled  or  if 
the  FSM  does  not  have  its  center  of  rotation  about  the  front  surface  of  the  mirror,  the  offset  at  the  focus¬ 
ing  lens  will  be  a  function  of  the  angular  disturbance.  For  example,  if  the  angular  disturbance  has  an 
incidence  angle  6.  with  respect  to  the  FSM  normal  when  it  is  at  null  and  the  mirror  normal  is  displaced 
by  an  angle  0^  from  null,  then  if  the  center  of  rotation  of  the  FSM  is  a  distance  D  behind  the  front 
surface  of  the  mirror  it  can  be  shown  that 


0  =  2Dsin  (6.  +e^)[l  -cosO^  )] 


(23) 


If  the  6^  is  the  nominal  angle  of  incidence  and  the  mirror  operates  around  its  null,  then 

O  =  2  D  Sin  Oo)  (24) 

Since  ©^represents  one-half  the  estimated  incidence  angle,  it  is  straightforward  to  show  that,  unless  the 
standard  tracking  algorithm  is  modified,  the  maximum  rejection  is  given  approximately  by 

R^,  =  Dsin(eo)(l-l)  .  (25) 

If  D  =  20  mm,  6^  =  45®,  f  =  20  cm,  and  L  is  positioned  to  within  1  percent,  then  R^^^x  • 


4.  EXPERIMENTAL  RESULTS 


The  goal  of  our  experimental  work  was  to  demonstrate  two-axis  high-precision  tracking  using  rel¬ 
atively  large  bandwidth  and  low  signal  power.  A  block  diagram  of  our  experimental  configuration  is 
depicted  in  Figure  6.  The  light  source  was  a  GaAlAs  laser  diode  operating  at  0.83|im.  Its  optical  output 
was  sinusoidally  intensity  modulated  at  54  kHz  with  a  modulation  depth  that  could  be  varied  from  0  to 
100  percent.  An  anamorphic  pair  was  used  to  convert  the  light  beam  from  elliptical  to  circular.  A  20: 1 
beam  expander  and  an  aperture  stop  were  used  to  transform  the  Gaussian  intensity  distribution  to  a 
5-mm-diam.  collimated  beam  of  uniform  intensity.  This  beam  was  directed  onto  a  beam-steering  mirror 
that  was  used  to  introduce  two-dimensional  angular  disturbances  into  the  tracking  loop.  Simulated 
platform  disturbances  as  well  as  sinusoidal  disturbances  could  be  injected  to  analyze  loop  performance. 
The  beam-steerer  had  internal  angular  monitoring  so  that  the  disturbance  injected  into  the  loop  could  be 
accurately  determined. 

A  1 : 1  relay  was  used  to  relay  the  pupil  from  the  randomizing  beam-steerer  onto  the  FSM.  If  the 
loop  was  operating  correctly,  the  FSM  would  correct  for  the  disturbance  introduced  by  the  randomizing 
beam-steerer.  The  light  was  then  split  into  three  different  paths.  One  path  led  to  a  television  monitor. 
The  focusing  lens  and  imaging  optics  were  configured  so  that  the  Airy-disk  pattern  could  be  viewed. 
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Figure  6.  Experimental  setup. 
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This  configuration  allowed  for  easy  alignment  and  a  coarse  check  of  tracking  performance.  The  middle 
path  served  as  a  reference  path  to  monitor  tracking  performance.  The  reference  error  signal  was  formed 
by  using  a  commercially  available  PIN  quadrant  detector,  four  transimpedance  amplifiers  (TlAs),  and  a 
summing  network.  The  summing  netw'ork  derived  the  baseband  azimuth  and  elevation  error  signals  by 
simple  addition,  subtraction,  and  low-pass  filtering  to  remove  the  IF  components. 

The  top  path  formed  part  of  the  tracking  loop.  A  neutral  density  filter  serv  ed  to  set  the  nominal 
received  signal  power.  To  achieve  low-noise  performance,  a  quadrant  APD  was  needed.  Unfortunately, 
there  are  no  commercially  available  quadrant  APDs  that  meet  our  requirements  of  low'  noise,  narrow' 
dead  zone,  and  gain  uniformity.  Therefore,  we  created  one  by  using  four  discrete  APDs,  some  imaging 
lenses,  and  a  four-way  reflecting  pyramid  splitter.  A  302-mm  focusing  lens  w'as  used  to  focus  a  spot 
onto  the  tip  of  the  image  splitter.  The  focal  length  was  chosen  to  form  a  spot  size  that  was  more  than  ten 
times  the  dead  zone  of  the  image  splitter.  Four  lenses  placed  one  focal  length  away  from  the  tip  of  the 
pyramid  were  used  to  collect  the  reflected  light  beams  and  image  a  pupil  onto  each  of  the  detectors  lo¬ 
cated  one  focal  length  away.  Imaging  a  pupil  onto  each  of  the  four  detectors  minimized  the  motion  of 
the  spot  across  the  detector  surface,  and  therefore  minimized  the  effect  of  gain  uniformities  across  the 
face  of  the  APD. 

Each  of  the  four  APDs  had  a  slightly  different  gain  dependence  on  bias  voltage.  This  variation, 
combined  wdth  the  open-loop  accuracy  of  the  bias  voltage  supplies,  led  to  the  need  for  a  gain-balancing 
algorithm.  Each  APD  would  be  biased  according  to  the  desired  gain  and  gain-v.s-bias-voltage  curve  for 
APD  1.  The  light  hitting  APD  1  was  maximized  by  monitoring  the  voltage  at  54  kHz  from  TIA  1  using 
a  hill-climbing  procedure.  Next,  the  light  out  of  APD  2  was  maximized.  The  bias-voltage  supply  for 
APD  2  was  then  servoed  until  the  voltage  out  of  TIA  2  was  equal  to  TIA  1.  The  remaining  APDs  were 
biased  in  the  same  fashion.  This  procedure  was  automated  on  a  computer  controller.  Successful  results 
were  obtained.  In  fact,  if  the  bias  voltages  corresponding  to  two  gain  settings  were  recorded  before 
tracking  was  initiated,  it  was  found  that  the  gains  could  be  switched  during  tracking  without  losing  lock. 
This  procedure  also  canceled  the  effects  of  gain  differences  in  some  of  the  front-end  electronics. 

The  outputs  of  the  APDs  were  fed  directly  into  low-noise  TIAs.  A  10-Mf2  feedback  resistor  was 
used  in  each  TIA.  This  value  represented  a  good  compromise  between  minimizing  thermal  noise,  the 
effect  of  amplifier  voltage  noise,  and  allowing  a  front-end  bandwidth  in  excess  of  100  kHz.  At  54  kHz, 
the  total  effective  front-end  noise  was  only  1  dB  above  the  thermal  noise  of  the  feedback  resistor.  The 
TIA  outputs  were  high-pass  filtered  to  remove  low-frequency  noise,  amplified,  and  sent  into  a  summing 
network  to  derive  azimuth,  elevation,  and  sum  signals.  The  three  signals  were  bandpass  filtered  and  then 
sent  into  a  normalizing  network.  The  bandwidth  of  the  sum-channel  filter  could  be  made  relatively 
narrow  since  it  does  not  affect  loop  dynamics  (except  during  changes  in  signal  power).  The  bandwidths 
of  the  difference  channels  must  be  wide  relative  to  the  desired  loop  bandwidth  so  that  the  loop  phase 
margin  is  not  compromised,  but  narrow  enough  so  that  the  noise  at  the  harmonics  of  the  carrier  is  elimi¬ 
nated  before  mixing  down  to  baseband. 

To  achieve  constant  loop  bandwidth  over  a  wide  range  of  received  signal  power,  it  was  necessary 
to  normalize  the  signals  proportional  to  the  angular  errors  by  the  power  in  the  sum  channel.  The  sum 
channel  was  split  into  two  paths.  The  first  path  was  sent  to  a  magnitude  detector,  the  output  of  which 
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was  used  in  a  iranslinear  circuit  to  normalize  the  difference  signal.  The  second  path  was  sent  to  a  hard 
limiter  that  was  used  to  create  a  reference  to  demodulate  the  normalized  difference  channels.  Normaliz¬ 
ing  at  IF  and  then  mixing  down  to  baseband  as  opposed  to  demodulating  and  then  normalizing  has  the 
advantage  of  minimizing  the  effect  of  dc  offsets  since  the  system  can  be  ac  coupled  right  up  to  the 
demodulator.  Identical  designs  are  used  for  the  azimuth  and  elevation  loops.  (In  addition  to  this  nor¬ 
malization  technique,  another  based  on  converting  the  angular  error  into  a  relative  phase  shift  between  a 
pair  of  IF  signals  and  then  using  double-balanced  phase  detectors  to  extract  the  normalized  angle  error 
was  investigated.^^  The  two  methods  yielded  similar  results.) 

The  normalized  error  signal  is  sent  into  the  loop  compensation  network.  The  compensator  w'as 
designed  to  achieve  a  type-two  loop  with  a  crossover  frequency  of  1  kHz  and  to  provide  6  dB  of  gain 
margin  and  45°  of  phase  margin  at  the  crossover  frequency.  The  open-loop  transfer  function  from  the 
input  to  the  power  amplifier  to  the  output  of  the  angle  error  detector  is  indicated  in  Figure  7.  This 
transfer  function  was  measured  by  driving  the  input  to  the  FSM  power  amplifier  and  measuring  the 
output  of  the  loop-angle  error  detector.  Note  that  the  FSM  behaves  as  a  lightly  damped  second-order 
circuit  with  a  natural  frequency  of  about  1 12  Hz.  This  frequency  is  equal  to  the  square  root  of  the  flexure 
spring  constant  divided  by  the  moment  of  inertia  of  the  rotor  and  mirror.  A  secondary  resonance  is 
evident  at  about  9  kHz.  It  is  this  resonance  that  limits  the  ultimate  closed-loop  bandwidth  obtainable 
with  the  FSM. 


Figure  7.  Measured  open-loop  transfer  functions. 
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The  departure  from  180‘'  of  the  uncompensated  transfer  function  is  mainly  a  result  of  losses  in  the 
magnetic  material  of  the  motor.  This  loss  acts  like  a  resistor  in  parallel  with  the  actuator  coils  and  gives 
rise  to  an  additional  15"  of  phase  shift  at  1  kHz.  This  negative  phase  shift  in  the  loop  transmission  is 
distinct  from  the  phase  shift  due  to  the  inductance  of  the  beam-steering  coil,  which  was  reduced  to 
approximately  1"  at  crossover  by  means  of  a  power  amplifier  with  extremely  high  output  impedance  to 
frequencies  in  excess  of  50  kHz.  The  design  combines  feedback  to  establish  output  current  level  w  ith  an 
open-loop  common-base  buffer  amplifier  to  provide  the  high  output  impedance. 

To  achieve  design  objectives,  the  compen.sator  included  two  integrators  to  increase  low-frequenc) 
loop  gain,  two  complex  zeros  located  at  1 10  Hz  with  a  damping  ratio  of  0.25,  a  zero  at  400  Hz  to  intro¬ 
duce  positive  phase  shift  at  crossover,  and  a  second-order  low-pass  filter  with  a  damping  ratio  of  0.707 
located  at  4  kHz.  This  low-pass  filter  reduced  noise  bandwidth  and  lowered  loop  gain  at  the  frequency  of 
the  secondary  resonance  and,  most  important,  added  more  than  90"  of  phase  shift  to  make  the  resonance 
phase  stable.  In  addition,  a  lead  network  furnishing  about  20"  of  positive  phase  shift  at  1  kHz  was  used 
to  compensate  the  negative  phase  shift  introduced  by  the  lossy  magnetic  material.  The  loop  transmission 
that  results  when  this  compensator  is  combined  with  the  fixed-element  dynamics  shown  earlier  is  also 
indicated  in  Figure  7.  This  figure  was  derived  from  the  measured  closed-loop  rejection  transfer  function. 

The  uncompensated  tracking  error  can  be  found  from  the  closed-loop  rejection  transfer  function. 
This  function  was  measured  using  three  different  methods:  (1)  inject  an  error  signal  from  inside  the  loop 
and  measure  the  resulting  tracking  error  internally;  (2)  inject  a  disturbance  from  outside  the  loop  using 
the  beam-steerer  and  measure  the  resulting  tracking  error  internally;  and  (3)  inject  a  disturbance  from 
outside  the  loop  using  the  beam-steerer  and  measure  the  resulting  tracking  error  externally  using  the 
reference  detector.  For  frequencies  above  1  Hz,  all  three  methods  yielded  the  same  results.  However, 
below  1  Hz  only  the  second  method  had  sufficient  resolution  to  obtain  useful  measurements  below 
-120  dB.  A  possible  reason  that  method  (3)  did  not  yield  better  low-frequency  rejection  measurements 
is  because  of  1/f  noise  in  the  reference  detector.  The  results  are  shown  in  Figure  8.  The  closed-loop 
transfer  function  is  also  indicated  in  this  figure.  The  angular  disturbance  of  the  beam-steerer  at  0. 1  Hz 
was  approximately  1  mrad  and,  therefore,  the  135  dB  of  rejection  corresponds  to  a  tracking  error  of  less 
than  1  nrad. 

The  NEA  as  a  function  of  signal  power  was  measured  by  normalizing  the  rms  voltage  out  of  the 
reference  detector  by  the  angular  gain  of  the  detector  with  the  disturbance  generator  turned  off.  The 
angular  gain  was  accurately  determined  using  two  different  methods.  Tht  first  method  involved  a  linear 
least-square  fit  to  the  measured  central  ±0.4  beamwidth  of  the  angle  discriminator.  The  second  method 
involved  inserting  a  precision  optical  wedge  into  the  beam  between  the  FSM  and  the  first  beamsplitter. 
The  two  methods  agreed  to  within  5  percent.  The  average  of  the  two  measurements  was  used  in  subse¬ 
quent  analyses. 

NEA  was  recorded  for  a  variety  of  signal  powers  and  APD  gains.  The  results  are  plotted  in  Fig¬ 
ure  9.  Our  NEA  data  were  measured  on  an  earlier  design  that  featured  a  500-Hz  crossover  frequency, 
and  the  vertical  axis  on  the  right  corresponds  to  those  measurements.  The  vertical  axis  on  the  left  is 
the  equivalent  NESD.  This  was  derived  from  the  measured  NEA  using  the  measured  noise  bandwidth 
(2200  Hz).  The  noise  bandwidth  was  found  by  numerical  integration  of  the  measured  closed-loop 
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Figure  8.  Measured  closed-loop  transfer  functions. 


Figure  9.  NESD  and  NEA  vs  APD  gain. 
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transfer  function.  TTie  symbols  represent  the  measured  azimuth  and  elevation  tracking  errors,  and  the 
lines  represent  theoretical  predictions  based  on  the  parameters  in  Table  1.  Sinusoidal  intensity  modula¬ 
tion  with  a  modulation  depth  of  0.5  was  used,  and  -1.5  dB  of  implementation  loss  was  assumed.  Note 
that  even  with  a  noise-induced  tracking  error  as  low  as  0.003  beamwidth,  the  equations  developed  earlier 
for  NESD  still  accurately  predict  the  performance. 

It  is  of  interest  to  compare  the  best  tracking  performance  achieved  at  each  of  the  three  signal 
powers  with  the  tracking  performance  bound  of  Equation  (10).  The  best  measured  NESD  values 
were  3.84  x  10'^,  4.16  x  10'^,  and  1.59  x  10'^  BW^/Hz  for  signal  powers  of  5.3  x  10^,  5.6  x  10^,  and 
4.2  X  10^  p/s,  respectively.  The  corresponding  losses  from  the  optimal  value  are  easily  calculated  to  be 
10.0,  10.6,  and  15.2  dB,  respectively.  Accounting  for  the  suboptimal  performance,  we  find  4.8  dB  due  to 
discriminator  gain  (Section  2.4.),  3  to  5  dB  due  to  device  noise  (Section  2.5.),  1 .5  dB  due  to  implementa¬ 
tion  loss  (above  paragraph),  and  1.0  dB  in  quantum  efficiency  (see  Table  I).  Note  that  an  additional  loss 
of  approximately  2.9  dB  for  the  4.2  x  10^  p/s  case  is  most  likely  accounted  for  by  the  fact  that  the  opti¬ 
mum  gain  was  above  100,  as  seen  in  Figure  8. 

The  system  was  designed  to  operate  over  a  static  and  dynamic  signal  power  range  of  100  to  1 . 
Various  tests  showed  that  the  system  performed  quite  well  over  this  range.  In  addition  to  the  static 
dynamic  range  capability  displayed  in  Figure  9,  we  investigated  the  tracking  performance  in  the  pres¬ 
ence  of  dynamic  fades.  We  found  that,  in  the  presence  of  smooth  or  instantaneous  fades,  the  NEA  would 
simply  increase  from  the  high-signal-power  value  to  the  low-signal-power  value.  We  also  measured 
tracking  system  bias  as  a  function  of  signal  power,  and  found  that  the  bias  was  limited  to  less  than 
0.05  beamwidth.  The  addition  of  a  simple  AGC  circuit  between  the  normalizing  circuit  and  the  sum- 
and-difference  channels  could  be  used  to  increase  the  dynamic  range  still  further. 

5.  CONCLUSIONS 

We  have  developed  a  variety  of  bounds  for  the  noise-equivalent  spectral  density.  Under  signal 
shot-noise-limited  operation,  direct  detection  was  shown  to  be  able  to  outperform  a  heterodyne  detection 
system  by  3  dB.  Also  under  signal  shot-noise-limited  operation,  the  performance  of  a  quadrant  detector 
was  within  4.8  dB  (2  dB)  of  being  the  optimal  detector  array  for  an  Airy-disk  (Gaussian)  intensity  dis¬ 
tribution.  As  the  background  and  device  noise  increase,  the  performance  difference  approaches  zero. 

The  optimal  quadrant  detector  radius  was  approximately  equal  to  the  first  null  in  the  Airy-disk  profile  for 
background-limited  operation.  If  the  tracking  detector  array  was  optimized  for  centroid  estimation,  near- 
optimal  tracking  performance  could  still  be  achieved.  Based  on  measured  parameters,  we  showed  that  an 
APD-based  tracking  system  can  outperform  a  PIN-based  tracking  system  by  as  much  as  20  dB. 

The  effects  of  non-focal-plane  processing  and  focus  errors  on  tracking  performance  were  quantified 
in  terms  of  the  effect  on  the  NESD,  system  rejection,  and  dynamic  range.  We  showed  that  to  keep  the 
loss  in  noise  performance  acceptably  small,  the  focus  error  should  be  kept  to  less  than  \/13  and  that  the 
detector  array  should  be  positioned  to  within  the  diffraction-limited  spot  size  times  the  system  f-number. 

Our  experimental  results  showed  very  good  agreement  with  theoretical  predictions.  Over  2.4  kHz 
of  closed-loop  bandwidth,  134  dB  of  rejection  at  low  frequencies,  and  tracking  errors  as  low  as  3/1000*^ 
of  a  beamwidth  were  demonstrated.  Low-noise  performance  was  demonstrated  using  four  single¬ 
element  APDs  and  an  image  splitter. 
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